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The influence of defect dipoles on the electrocaloric effect (ECE) in acceptor doped BaTiOs 
is studied by means of lattice-based Monte-Carlo simulations. A Ginzburg-Landau type effective 
Hamiltonian is used. Oxygen vacancy-acceptor associates are described by fixed defect dipoles 
with orientation parallel or anti-parallel to the external field. By a combination of canonical and 
microcanoncial simulations the ECE is directly evaluated. Our results show that in the case of 
anti-parallel defect dipoles the ECE can be positive or negative depending on the density of defect 
dipoles. Moreover, a transition from a negative to positive ECE can be observed from a certain 
density of anti-parallel dipoles on when the external field increases. These transitions are due to 
the delicate interplay of internal and external fields, and are explained by the domain structure 
evolution and related field-induced entropy changes. The results are compared to those obtained 
by MD simulations employing an ab initio based effective Hamiltonian, and a good qualitative 
agreement is found. In addition, a novel electrocaloric cycle, which makes use of the negative ECE 
and defect dipoles, is proposed to enhance the cooling effect. 

PACS numbers: 77.70.-fa, 77.80.-e, 77.80.Jk, 74.62.Dh 


I. INTRODUCTION 

Ferroelectric materials based on BaTiOs are of inter¬ 
est for solid state refrigeration because of their significant 
electrocaloric effect (ECE).^^® Moya et al.^ measured the 
ECE of single crystalline BaTiOs samples in the direct 
and indirect way and reported temperature variations 
of around 1.0 K. Qian et al.^ investigated the ECE in 
BaZro, 2 Tio. 80 s and found a large temperature variation 
of 4.5 K over an operating temperature range of 30 K, 
which is wider than that in BaTiOs. Recently, experi- 
mentaR^® and theoreticaR*^’^^ studies revealed the coex¬ 
istence of the positive (AT > 0) and negative (AT < 0) 
ECE. They concluded that if the external field E®’' is not 
parallel with the dielectric polarization of the material, 
may increase the amount of dipolar entropy and thus 
reduce the phonon contribution under isentropic condi¬ 
tions causing a negative AT. In first-principles based 
calculation Ponomareva and Lisenkov^° found that if the 
direction and strength of the external field are varied the 
sign of the ECE can be efficiently controlled near the 
phase transition temperature. Pirc et al}^ studied the 
negative ECE for an anti-ferroelectric materials based on 
a generic Kittel model and found transitions from a nega¬ 
tive to a positive ECE depending on E®^ below the Curie 
temperature. Axelsson et al}^ modeled the ECE by di¬ 
rectly evaluating the entropy from the partition function 
of a one-dimensional lattice model and showed how the 
field applied along a particular lattice direction induces 
a negative ECE. According to their model the sign re¬ 
versal of the ECE happens when two phase transition 
temperatures are close. 

Even though the influence of phase transitions on 
the ECE in ferroelectrics and ferroelectric relaxors was 


widely studied, there is very few work to reveal the role 
of dopants and defects on ECE. 

In acceptor doped BaTiOs, Ba or Ti ions are substi¬ 
tuted by ions with a lower valence. In this case, charge 
neutrality is typically obtained by compensating oxygen 
vacancies. Acceptor doped materials, which are also re¬ 
ferred to as hard doped, are difficult to polarize, exhibit 
high coercive fields as well as small strains.It was well 
established that in acceptor doped materials defect reac¬ 
tions between dopants and oxygen vacancies can occur. 
Eor BaTiOs it was demonstrated that (Mn^(.j — as¬ 

sociates with an excess orientation parallel to the spon¬ 
taneous polarization are present in reduced crystals.^® 
Non-switchable defect dipoles impose a restoring force 
for reversible domain switching^® and thus might have an 
impact on the ECE, since their gradual re-orientation is 
determined by the barrier for oxygen vacancy migration, 
hence these defect complexes cannot immediately follow 
the polarization switching. The possibility to enhance 
and control the ECE in dielectric materials by the pres¬ 
ence of internal dipoles was formulated by Van Vechten^® 
in an US-patent in the late seventies. Only recently, how¬ 
ever, Griinebohm et al}^ showed by means of Molecular 
Dynamics (MD) simulations based on an effective Hamil¬ 
tonian approach that the ECE can switch in the presence 
of fixed defect dipoles from positive to negative behavior. 

In this work, we simulated the influence of the defect 
dipoles on the ECE, using the MC model presented in our 
previous work^°. For the sake of completeness the model 
is shortly introduced in Sec. II. The results are presented 
and discussed in Sec. Ill with three different subsections. 
In Subsec. Ill A the influences of different types of de¬ 
fect dipoles on the ECE are investigated. Results show 
that the ECE at high temperature can be increased in 
presence of defect dipoles which lie parallel to the exter- 
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nal field (parallel defect dipoles). Moreover, a negative 
ECE and a double-peak behavior can be induced by de¬ 
fect dipoles which are anti-parallel to the external field 
(anti-parallel defect dipoles). In Subsec. IIIB we com¬ 
pare the ECE results by MC and MD simulations for 
different loading scenarios (field-on / field-off). It shows 
good qualitative agreements on the peculiar ECE behav¬ 
iors, even though in the MC simulations the strain cou¬ 
pling is absent. In Subsec. IIIC strategies to tailor the 
ECE cycle by defect engineering are discussed. By apply¬ 
ing additionally a reversed electric field during cooling, 
an enhanced temperature drop is observed in the sample 
with 3% defect dipoles. Finally, the results are summa¬ 
rized in Sec. IV. 


II. MODEL AND SIMULATION SETUP 

In this study, we analyze the impact of fixed parallel 
and anti-parallel defect dipoles on the ECE by directly 
simulating the adiabatic temperature changes using 2 - 
D lattice-based MC simulations based on a Ginzburg- 
Landau type Hamiltonian for BaTiOa of our previous 
work (see Ref. 20). The parameters used in this paper are 
identical to the ones given in Ref. 20. For clarification, 
the model is shortly summarized here. The total energy 
of the system E* involves the thermal energy Ek and the 
potential energy E: 

E*=Ei, + E. (I) 

The potential energy E is assumed to be of a Ginzburg- 
Landau type, which includes four contributions: the Lan¬ 
dau multi-well energy term Ed, the dipole-dipole interac¬ 
tion energy E^ip, the domain wall energy Eg,, that arises 
from short-range and elastic interactions,^^ and the elec¬ 
trostatic energy Eg: 

E = Ed + Edip -|- Egr -|- Eg. (2) 

In contrast to the previous work^°, instead of random 
fields, defect dipoles are considered. The model is ex¬ 
tended by introducing fixed non-switchable polarizations 
on certain lattice sites to mimic the presence of defect 
dipoles. The elastic energy and the electrostrictive inter¬ 
action are not explicitly included in the Hamiltonian^^, 
but the elastic energy contribution to domain wall en¬ 
ergies is implicitly considered in the gradient term Eg,. 
However, this simplification has no qualitative influence 
on the calculated entropy variations, since defect dipoles 
have only a weak elastic strain field and the possible im¬ 
pact of elastic energies on the domain switching behavior 
is the same in the isothermal and adiabatic configurations 
with applied field. This statement has been confirmed by 
a comparison with MD simulations taking the local strain 
into account (see Subsec. IIIB). 

For the evaluation of the ECE by the MC simula¬ 
tions, there are two methods: the indirect method uti¬ 
lizing Maxwell relations, and the direct method^® using 
the modified Creutz’s algorithm. The accuracy of the 


indirect method depends significantly on the numerical 
integration.^^ By contrast, the direct method can sample 
the change of the temperature accurately, without de¬ 
pendence on the numerical integration.^® In this paper 
the direct method is adopted. 

The reason for aging in poled materials is known to 
origin from the alignment of the dipole defects^^ and is 
thus important for possible applications. However, since 
the orientation of defect dipoles follows the spontaneous 
polarization on time scales of days to months depend¬ 
ing on temperature, and the simplifying assumption of 
a fixed dipole orientation within an ECE cycle occur¬ 
ring on much shorter time scales is reasonable. Unidi¬ 
rectional defect dipoles might be introduced through a 
long-time poling process.Defect dipoles with differ¬ 
ent orientations are investigated in this paper, namely 
defect dipoles pointing parallel to E®’' and defect dipoles 
oriented anti-parallel to E®’'. The local polarization of 
the defect dipoles^® can be approximated as Ed = Q^/Vo, 
where the charge of the neutral defect dipole in BaTiOa 
is <7 = -|- 2 e, I is half of the lattice constant oq and Vq is 
the lattice volume. Therefore, the local polarization for 
the defects is assumed as I.0Cm“^, and the direction of 
these defect dipoles is fixed. When several lattice sites are 
occupied by one type of defect dipole, i.e. the parallel or 
anti-parallel defect dipoles, the inhomogeneous internal 
field Ej becomes stronger, respectively. Thus, an impact 
on the accessible entropy changes can be expected. 

It should be noted rather high fields up to 
122.2kVmm“^ are applied in this model. However, ex¬ 
trinsic effects, which is caused by localized nucleation of 
domains with reversed polarization in real ferroelectrics, 
can significantly reduce the coercive field in larger sam¬ 
ples or polycrystalline materials. 

Samples with high density of the defect dipoles 
can be synthesized, e.g. Ba(Tio. 93 Feo. 07 ) 03 ®^ and 
Ba(Tio. 95 Mno. 05 ) 03 ®^ thin films, and the defect density 
studied in our model (< 6 %) is within this justifiable 
region. 


III. RESULTS AND DISCUSSION 
A. Positive and negative electrocaloric effect 

For the simulations in this subsection the following 
loading history is applied. Firstly, the sample is equi¬ 
librated at a given temperature without applying a ex¬ 
ternal field for 28,800 MC steps. During one MC step all 
L X L sites are visited, and the polarization switching is 
allowed with a certain probability. From this equilibrium 
state the sample is then poled at a fixed external electric 
field for 28,800 MC steps. After instantaneous removal 
of the poling electric field and a third relaxation step 
with 28,800 MC steps, the prepoled sample is obtained. 
These processes are performed in a canonical ensemble 
at constant temperature Tq. Finally, the prepoled sam¬ 
ple is used to study the ECE under an adiabatic condition 
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FIG. 1. Positive and negative ECE. (a) ECE in the presence of parallel aligned defect dipoles. Different defect concentrations 
are considered under a given external field (b) ECE as a function of E®’' for a defect concentration of 6%. (c) Positive 

and negative ECE in presence of anti-parallel dipoles. When the defect density exceeds a critical value, the resultant internal 
anti-parallel field overcomes E®’‘ and a negative temperature change is observed (negative ECE). (d) Influence of E®’®, while 
anti-parallel defect dipoles with a concentration of 3% are present. When E®^ surpasses the internal field Ei induced by the 
anti-parallel defect dipoles, a negative ECE appears only at lower temperatures, while a positive ECE dominates at higher 
temperatures, see the olive line (84.6kVmm“^). The above phenomena are explained by the domain structures at the points 
I, II, III, A and B in Fig. 2. 


with a field switched on instantaneously. Hence the mi- 
crocanonical ensemble with a constant total energy is uti¬ 
lized. After running 43,200 MC steps, the system reaches 
the equilibrium state with temperature Ti, i.e. the end 
of the adiabatic stage. This kind of loading is named the 
“Field on” case. The typical measurement time length of 
the ECE in experiments is seconds, which is much longer 
than the characteristic switching time length of the po¬ 
larization. Therefore, the instantaneous field removal or 
field application is utilized. Although the ramping of the 


external field can give different values for the temperature 
change, it has been shown that the resultant difference 
is small.In Fig. I the ECE is investigated within the 
“Field on” case. 

The influence of the concentration of parallel defect 
dipoles on the ECE is illustrated in Fig. 1(a). For 
pure BaTiOa the magnitude of AT/AE at the peak is 
5.7 X 10“^Km/V, which is comparable with the exper¬ 
imental data 8.3 x 10“^Km/V by Moya et al?. With 
increasing defect concentration, the ECE peak decreases. 
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FIG. 2. Domain structure explanation of ECE for the points I, II, III, A and B marked in Fig. 1. The external field is 
applied in the horizontal direction and points to the right, (a) and (b) are the domain structure snapshots before and after the 
adiabatic stage for the point I; (c) and (d) for point II; (e) and (f) for point III; (g) and (h) for point A; (i) and (j) for point 
B. The black dots with white arrows represent the defect dipoles, either parallel or anti-parallel to The red, blue, green 
and yellow dots represent the dipoles pointing respectively to the left, to the right, to the bottom and to the top. 


and the maximum shifts to higher temperatures, due to 
the reason that Ei induced by the defect dipoles stabilizes 
the tetragonal phase. However, the interesting feature 
shown in Fig. 1(a) is that the temperature change of the 
samples with parallel defect dipoles is larger than that 
in the defect-free one within the high temperature range. 
In addition, the obtainable ECE is enhanced when E®’' 
is increased (see Fig. 1(b)) since higher E®’' can induce a 
more ordered system. Similar results have been reported 
by Rose and Cohen^^, where no defect dipoles were con¬ 
sidered. It can be concluded that above the phase tran¬ 
sition temperature the ECE can be enhanced either by 


application of higher fields or by incorporation of parallel 
defect dipoles. 

For clarification of the observed effects we study the 
domain configurations at different points in the cycle. 
The prepoled sample shows a multiple domain configu¬ 
ration (Fig. 2(a)) which is switched into a state of reduced 
dipolar entropy after the held is switched on (Fig. 2(b)). 
And thus the temperature is increased under adiabatic 
conditions by E®’'. It can clearly be seen that a smaller 
number of domains occurs in the presence of E®’' than in 
the prepoled sample without E®’'. 

If we compare now with the case, where defect dipoles 
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in parallel orientation to E®’' are present, the situation 
changes (Fig. 2(c) and (d)). In the initial state the con¬ 
figurational space is reduced as compared to the defect 
free sample (see Fig. 2(a) and (c)), since defect dipoles 
induce E, that locally stabilizes the domain structure 
against thermal fluctuations. If E®’' is now switched on, 
the change in configurational entropy (^conf) is smaller 
than in the defect-free case and thus a smaller temper¬ 
ature variation is possible, which explains the defect- 
concentration dependence of the ECE. If fewer domain 
walls disappear, the decrease of the domain wall energy 
is reduced and the excess heat is reduced. In contrast to 
the low temperature range, at higher temperatures the 
ECE is elevated by the existence of the parallel defects. 

A different scenario occurs, if the defect dipoles are 
pre-aligned in anti-parallel direction (see Fig. 1(c)); a 
situation which can be installed by poling the sample on 
extended time scales, before reversing the field direction. 
For defect concentrations above 2% the ECE turns from 
a positive effect into a negative one and becomes more 
pronounced with increasing E®’'. This can be explained 
again from the characteristic domain structure shown in 
Fig. 2(e) and (f) with 3% anti-parallel defect dipoles. Ini¬ 
tially, the system is close to a single domain configuration 
(Fig. 2(e)) with orientation parallel to the defect dipoles 
and thus has a small S'conf because of the large average 
internal field of (Ei) = —87.2kVmm“^. 

If now an external field of E®’^ = 47.0kVmm“^ is 
switched on, Ei is partly compensated and the mate¬ 
rial can locally sample a multi-domain structure with in¬ 
creased S'conf (see Fig. 2(f)) corresponding to a reduced 
potential energy. Thus, vibrational entropy has to be re¬ 
moved from the phonons under isentropic conditions and 
the negative ECE appears. It should be noted that the 
domain snapshots, i.e. Fig. 2(e) and (f), are for 400 K. At 
the negative ECE peak, i.e. 360 K, the prepoled sample 
is more like a single domain type, and the ^conf change 
is more prominent. 

When the magnitude of E®’' is, however, compara¬ 
ble with that of E, (hereby between 84.6kVmm“^ and 
112.8kVmm“^), both the ECE and the negative ECE 
start to coexist (see Fig. 1(d)). A sharp rise in tem¬ 
perature is observed at the transition from the negative 
into positive temperature change, and both the negative 
ECE peak and the first positive ECE peak shift to lower 
temperature with increasing E®’'. 

This transition, again, can be explained by inspect¬ 
ing the domain structure for the case of E®’^ = 

112.8kVmm“^. Before E®’' is applied, thermal energy 
fluctuations provide the ability to sample configurational 
space in areas, which are not pinned by the presence 
of local defect dipoles. This can be seen in Fig. 2(g), 
where the single domain states (red) form a percolat¬ 
ing (non switching) network of the sites localized in the 
direct vicinity of defect dipoles, while multiple domain 
configurations, which can take various iso-energetic con¬ 
figurations, occur in between. If now an E®’', which is 
stronger than the internal field, is switched on, the situa¬ 


tion is reversing. Sites being part of the initially inactive 
percolating network around the anti-parallel defect sites 
(see Fig. 2(h)) are those where the external field is com¬ 
pensated. Thus, only these sites in the direct vicinity of 
defect associates can effectively sample configurational 
space, while the remaining matrix polarizations orient 
parallel to E®^ and do not contribute to S'conf- Since 
the relative number of sites which can possibly flip their 
configuration scales with the defect concentration and is 
thus smaller than the number of sites dominated by E®’', 
Sconf is reduced and the ECE becomes positive again. 

Under high enough field (hereby 103.4 and 
112.8kVmm“^) and within a low temperature range 
(< 240 K) the negative ECE still persists due to the 
same reason discussed above for the case under low 
fields. When the temperature slightly increases, a sharp 
transition from this negative to positive ECE is observed 
as well as in the case with 84.6kVmm“^. However, 
with further increasing of the initial temperature, 
AT/E in Fig. 1(d) shows a cusp (around 320 K under 
112.8kVmm“^) and a hump (around 400 K under 
112.8kVmm“^). This double-peak behavior can be 
attributed to a temperature induced phase transition in 
the initial state. As can be seen from the domain con¬ 
figuration, the local domains have formed at T = 400 K 
and the system has locally taken a multi-domain config¬ 
uration in the initial prepoled sample, which is higher in 
potential energy than the single domain case prevailing 
below T = 320 K. Thus the effective AT increases again. 
In brevity, with increasing initial temperature this local 
order is further reduced and the material behaves like a 
regular ferroelectric material. 


B. Comparison between MC and MD 

The coupling to the strain is not considered directly in 
the MC simulations. In order to justify this simplifica¬ 
tion and the validity of the results obtained, MD simu¬ 
lations in 3D have been performed deploying the feram 
code^^ based on an effective Hamiltonian. This Hamil¬ 
tonian includes the self energy due to the ferroelectric 
soft mode, the long-range dipole-dipole interaction en¬ 
ergy, the elastic energy, and the energy due to the cou¬ 
pling between the local soft modes and the strain. Details 
of this Hamiltonian were explained in Refs. 21, 34, and 
35. It should be noted that only three degrees of free¬ 
dom are explicitly taken into account in the MD simula¬ 
tions (the strain is optimized for each step but not used 
as degree of freedom). Thus, the obtained temperature 
change overestimates the adiabatic temperature change. 
The leading error can be corrected by a rescaling of the 
degrees of freedom with 3/15 (actual number of degrees 
of freedom/degrees of freedom of BaTiOa). The set of pa¬ 
rameters for the effective Hamiltonian for BaTiOa have 
been obtained by density functional theory simulation at 
T = 0 K and are listed in Ref. 36. 

Within the MD simulations, periodic boundary condi- 
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FIG. 3. Comparison of AT/E over T from the MC simulations with 3% anti-parallel defect dipoles (top row) with the results 
of MD with 1% anti-parallel defect dipoles (bottom row), (a) and (b) are “Field on” cases, while (c), (d), (e) and (f) are “Field 
off’ cases. For the “Field on” case, the ECE double-peak behavior is present in MC and MD simultaneously. For the “Field 
off’ case both in MC and MD, a single peak is observed under a low field ((c) and (d)), and a shape with one hump and one 
cusp is demonstrated under a high field ((e) and (f)). Therefore, a qualitative agreement between the MC and MD simulations 
exists. 


tions and a cell size of 38.4 nmx 38.4 nmx38.4 nm have 
been used. For simulations in the canonical ensemble 
(NPT), the Nose-Poincare thermostat is applied.De¬ 
fects are modeled by freezing the local soft mode on ran¬ 
dom positions to a certain value. In the following 1% 
polar defects anti-parallel to the external field are frozen 
to mimic a soft-mode amplitude of 0.2 A (approximately 
0.5 C m“^). 

The mean distance between defect dipoles in the 2D 
MC simulations is larger than that in the 3D MD for 
the same concentration. Namely, in 3D more neighbor¬ 
ing defect dipoles appear in the vicinity of non-defect 
dipoles than in 2D, which leads to a stronger influence of 
the defects. Hence, the defect concentration, as well as 
the strength of the defect dipoles, should be set smaller 
in the MD case to achieve similar effects of defects. In 
addition, the different energy contributions used in MC 
and MD, as well as the different parameterizations, re¬ 
sult in a quantitative difference in the coupling between 
the external field and the dipoles. Due to all these rea¬ 
sons, the strength of the electric field has to be differ¬ 
ent as well. Nevertheless, qualitative agreements can be 
reached, and even the quantitative agreement of AT/E 


is observed (see following). 

For comparison, MC and MD simulations are carried 
out for three different cases. The first is a “Field on” 
case with a strong field. The loading history of a ’’Field 
on” case has been explained in Subsec. Ill A. The other 
two are “Field off” cases with a weak and strong field, re¬ 
spectively. The so-called “Field off’ case involves the fol¬ 
lowing loading history. Firstly, the system is equilibrated 
at constant temperature under a fixed external electric 
field that is anti-parallel to the defect dipoles. Secondly, 
in MC the field is removed instantaneously, while in MD 
the field is ramped down with a rate of 0.001 kV/cm/fs. 
Simultaneously, the microcanonical ensemble is utilized 
at a constant energy, using the Multi-demon method in 
MC^° and the leapfrog method in MD. Hereby, the final 
state at the end of the stage with constant-temperature is 
used as the initial state for the microcanoncial ensemble. 
Thirdly, after the system reaches equilibration, the ther¬ 
mal energy is obtained so that the temperature change 
can be evaluated. In MD the kinetic energy is monitored 
after 500,000 further equilibration steps. By this slow 
ramping and long equilibration the system can overcome 
potential energy barriers due to domain switching and 
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is not stuck in a meta-stable state. Finally the temper¬ 
ature can be calculated by sampling the kinetic energy 
spatially and historically in both and MD simula¬ 

tions. A time step of At = 1 fs in MD is used in both 
ensembles. It is worth mentioning again that the results 
produced by the slow ramping are only slightly different 
from those by the instantaneous ramping.^® Therefore, 
the results are comparable between MC and MD even 
though different loading histories are utilized. 

It should be noted that the phase transition from 
the paraelectric to ferroelectric state in MC lies around 
400 K. However, in the MD simulations the ferroelectric- 
paraelectric phase transition temperature is systemat¬ 
ically underestimated (300 K) compared to the 400 K 
found in experiment.Nevertheless, the aim is to com¬ 
pare the results of MC and MD qualitatively rather than 
quantitatively. 

For the “Field on” case under a high field, both the MC 
and MD simulations demonstrate a double-peak ECE be¬ 
havior, as it is seen in Fig. 3(a) and (b). The reason for 
the double-peak has been explained in the discussion of 
Fig. 1(d). Note that in the MD simulations, the sample is 
not prepoled. Since the internal fields induced by the de¬ 
fect dipoles are fairly strong, after 500,000 equilibration 
steps the domain configurations before and after poling 
have small differences. 

For the “Field off” case under a weak field, both the 
ECE curves obtained by MC and MD simulations show 
a single peak as shown in Fig. 3(c) and (d). In the MC 
results the peak appears at 56.4kVmm“^, and in MD 
at 8.5kVmm“^. The underlying physics for this peak 
behavior has been explained in the discussion of Fig. 1(c) 
and (d). 

For the “Field off” case under a high field, the curve 
experiences a hump and a cusp, which is captured both 
in the MC and MD simulations (see Fig. 3(e) and (f)). 
The explanation of this peculiar ECE behavior can be 
referred again to the entropy change in this process. Be¬ 
low the phase transition temperature, the defect dipoles 
dominate the domain patterns, and the configuration is 
a single domain when no external fields is applied. At 
very low temperature (e.g. 120K), the external electric 
field induces disorder. After removing the fields the con¬ 
figurational entropy thus decreases slightly, which results 
in a small increase of the vibrational entropy, i.e. of the 
temperature. With increasing the initial temperature, 
i.e. the thermal fluctuations, it becomes much easier to 
trigger disorder under electric fields. Thus the change of 
the domain patterns after removing the external fields are 
more prominent, i.e. the decrease of the configurational 
entropy is higher during the removal of the fields. Hence, 
the increase of the vibrational entropy and of the temper¬ 
ature is more significant. This phenomenon corresponds 
to the temperature below 240 K in MC and below 180 K 
in MD. When the initial temperature is elevated further, 
due to the higher thermal fluctuations the domain pat¬ 
terns are more disordered in the sample both with and 
without fields. It signifies that the difference between the 


configurational entropy with and without fields becomes 
smaller. Therefore, the positive ECE weakens when the 
initial temperature is increased from 240 K to 400 K in 
MC and from 180 K to 280 K in MD. In other words, the 
positive ECE is largest at 240 K in MC and at 180 K in 
MD. When the temperature is high enough, the thermal 
fluctuation makes it fairly easy to reach an ordered state 
under assistance of the applied field. Then the disor¬ 
der increases if the external field is removed, which leads 
to an increase of the configurational entropy and a de¬ 
crease of the temperature. This phenomenon is common 
in both the MC and MD simulations above the corre¬ 
sponding para- to ferroelectric transition temperature. 
The upcoming cusp appears in MC at 440 K and in MD 
at 400 K. The existence of the cusp is due to the fact that 
the phase transition point is shifted to higher tempera¬ 
ture in the “Field off” case. This is similar to the shift of 
the transition temperature in presence of parallel defect 
dipoles for the “Field on” case in Fig. 1(a). 

In summary, the MC and MD simulations show good 
qualitative agreements in the three different loading 
cases. 


C. Defect engineering for the ECE 

Some researchers^°d3 have revealed a negative ECE in¬ 
duced by noncolinear-field-driven phase transition, and 
pointed out that this phenomenon can be utilized to en¬ 
hance the ECE. Similar to that concept, a parallel or 
anti-parallel electric field is applied with respect to the 
direction of the defect dipoles. Nevertheless, in this work 
the main focus of interest is on the influence of defect 
dipoles, as shown in the previous subsections. Based on 
the gained knowledge the strategies to tailor the ECE 
by defect engineering can be considered. It is of special 
importance since defects are unavoidable in real material 
systems. 

Fig. 1(a) indicates that at high temperature the sam¬ 
ples with defect dipoles can generate a higher ECE than 
that without defects under a parallel external field. For 
example, under an electric field of 47.0kVmm“^, the 
temperature change at 460 K for the sample with 2% par¬ 
allel defect dipoles is 12.1 K while only 7.2 K for the sam¬ 
ple without defects. In other words, the parallel defect 
dipoles have a potential for high temperature applica¬ 
tions. 

More importantly, the design of a novel electrocaloric 
cycle can be stimulated by utilizing the positive-negative 
ECE transition. As an example, we propose in the follow¬ 
ing a modified Carnot-like Cycle with a reversed external 
field, which provides the possibility to improve the cool¬ 
ing efficiency of the electrocaloric device. A cycle with 
3% defect dipoles at an initial temperature Tq = 400 K is 
illustrated in Fig. 4. In the classical concept the ECE cy¬ 
cle^® includes four states “a-b-c-d-a” (see dashed arrows), 
while the new concept contains four states ’a-b-c-e-a’ (see 
solid arrows). Compared with the state (d), the new state 
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FIG. 4. Proposed modified ECE cycle to enhance the resultant AT of a sample containing defect dipoles. The conventional 
cycle involves a-b-c-d-a (dashed arrows), and the new proposed one includes a-b-c-e-a (solid arrows). Under adiabatic condition 
a field is applied rapidly parallel to the defect dipoles, taking the initial state (a) at To to a state (b) with lower configurational 
entropy at Ti. By transferring heat to a heat sink, the state (b) transforms to the state (c) with further lower configurational 
entropy at Tq. Removing the field rapidly under adiabatic condition, from state (c) the system is cooled to state (d) at T 2 . By 
comparison, applying a field anti-parallel to the defect dipoles the system is cooling further from state (c) to state (e) (Ta < Ta). 
It should be noted that the magnitude of the external field should not be greater than the magnitude of the average internal 
field induced by the defect dipoles. Finally, by absorbing heat from a load, the system warms up to the state (a). The relation 
of the involved temperatures can be summarized as To < T 2 < To < T\. It can seen that through introduction of the defect 
dipoles, an electrocaloric cooling enhancement is achieved. 


(e) involves a reversed electric field, which should induce 
a further temperature drop, due to the negative ECE 
indicated in Fig. 1(d). 

The red and blue arrows represent the heating and 
cooling processes, respectively. Process i is an adiabatic 
process from state (a) with Tq = 400 K to (b) with 
Ti = 421.6 K, which decreases the configurational en¬ 
tropy, and increases the temperature by application of 
fields El. Hereby Ei is parallel to the defect dipoles. 
Process ii is an isothermal process from stage (b) to (c) 
with To = 400 K, in which the temperature decreases by 
releasing the heat to a heat sink. Process iii(l) from 
state (c) to (d) with T 2 = 380.9 K and iii(2) from state 


(c) to (e) with Ta = 371.5 K are both adiabatic processes. 
In process iii(l) the configurational entropy increases by 
removal of a field while in process iii(2) by application 
of a reversed field E 2 . Namely, the direction of E 2 is 
opposite to El. By application of this reversed field the 
vibrational entropy decreases more than in process iii(l). 
The temperature drop in iii(2) is 28.5 K, which is higher 
than 19.1 K in iii(l). Process iv(l) from state (d) to 
(a) and iv(2) from state (e) to (a) are both isothermal 
processes, in which the temperature increases though ad¬ 
sorption of heat from a heat source. The example shows 
that the incorporation of defect dipoles and application of 
a reversed field can enhance the cooling effect by around 
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50%. 


IV. CONCLUSION 

In summary, our MC simulation results reveal that the 
ECE can be influenced heavily by the presence of de¬ 
fect dipoles due to acceptor-vacancy associates, which 
agrees qualitatively with MD results. The incorporation 
of the parallel defects provides a application potential for 
high temperature by increasing the achievable tempera¬ 
ture change. The anti-parallel dipoles act as “memory 
elements” which reduce the configurational space in the 
prepoled sample and thus allow to turn the positive ECE 
into a negative ECE. Under certain fields both effects 
can be combined, which offers possibilities to increase 
the cooling efficiencies of electrocaloric devices, and indi¬ 


cates potential applications of ferroelectric materials with 
defect dipoles for the ECE. 
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